Effect of random interactions in spin baths on decoherence 
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We study the decoherence of a central spin 1/2 induced by a spin bath with intrabath interactions. 
Since we are interested in the cumulative effect of interaction and disorder, we study baths comprising 
Ising spins with random ferro- and antiferromagnetic interactions between the spins. Using the 
resolvent operator method which goes beyond the standard Born-Markov master equation approach, 
we show that, in the weak coupling regime, the decoherence of the central spin at all times is 
entirely determined by the local-field distribution or equivalently, the dynamical structure factor of 
the Ising bath. We present analytic results for the Ising spin chain bath at arbitrary temperature 
for different distributions of the intrabath interaction strengths. We find clear evidence of non- 
Markovian behavior in the low temperature regime. We also consider baths described by Ising 
models on higher-dimensional lattices. We find that interactions lead to a significant reduction of 
the decoherence. An important feature of interacting spinbaths is the saturation of the asymptotic 
Markovian decay rate at high temperatures, as opposed to the conventional Ohmic boson bath. 



PACS numbers: 

I. INTRODUCTION. 

Recent developments in nanophysics have made pos- 
sible the use of the charge and spin dynamics of elec- 
trons to develop new technologies like spintronics and 
quantronics. This has moreover led to the possibility of 
using the electron spin, or other more complex entities 
like the phase in Josephon junctions to fabricate qubits 
for quantum computing. However, the utility of these 
nanosystems as qubits is strongly limited by their cou- 
pling to the omnipresent dissipative environment. The 
environment destroys the coherence of the qubit over a 
certain time scale and a lot of recent theoretical and ex- 
perimental activity has focused on ways and means to 
increase this time scale. This clearly emphasizes the im- 
portance of understanding the effects engendered by the 
coupling of a two level system to a dissipative bath. 

The fact that the environment plays a crucial role 
in the physics of small quantum systems has been well 
known since the pioneering work of RefQ], where it was 
shown that the coupling of a two level system to an 
Ohmic boson bath could effectively suppress the tunnel- 
ing of the two level system. In the context of decoherence, 
the most commonly studied problem is the spin-boson 
models which describes the effect of a dissipative bosonic 
bath on a central spin, where the spin can be an effec- 
tive description of a system whose discrete lowest energy 
levels dominate the physics at low enough temperatures 
and the bosons are often the phonons present in the sys- 
tem. A physical manifestation of the spin boson problem 
is a nanomagnet (described by a giant spin) coupled to 
phonons 3 . However, for many practical realisations of a 
central spin or a qubit (spin 1/2), a spin bath comprising 
other spins might be the principal source of decoherence. 
This is indeed the case in semiconducting quantum dots, 
where the nuclei with non zero spins constitute the spin 
bath and interact with the central electronic spin in the 



dot via the hyperfine interaction^. Another manifesta- 
tion of a spinbath occurs in Si:P 6 . The abundance of spin 
baths in real systems, necessitates an understanding of 
their effect on decoherence. Unlike the case of bosonic 
baths often modeled as a collection of harmonic oscilla- 
tors, spin baths can exhibit a wide range of phenomena 
depending on the interactions between the spins, residual 
anisotropies etc. Clearly one expects any resulting deco- 
herence of the central spin to depend rather crucially on 
the underlying nature of the spin bath and its coupling 
to the central spin. 

Earlier studies which considered independent spins in 
the bath seemed to indicate that spin baths were not 
qualitatively different from bosonic baths^. More inter- 
estingly, recent studies of decoherence induced by spin 
baths described by mean field Hamiltonians have demon- 
strated that interactions between the bath spins can be 
used as a lever to augment the time scales over which 
the system decoheres2*i2*ii. These results were how- 
ever obtained either numerically for a bath with a small 
number of spins or for the special case where the bath 
Hamiltonian commutes with the bath-central spin cou- 
pling Hamiltonian leading to an effective classical deco- 
herence. A more robust treatment of intrabath interac- 
tions was presented in Refill where the authors studied 
numerically the zero temperature quantum decoherence 
of two coupled spins engendered by a bath described by 
the random transverse Ising model. The authors used 
this model to argue that the central spin decoheres dif- 
ferently, depending on whether the spinbath has a regu- 
lar or chaotic spectra. Despite their various drawbacks, 
these works collectively highlight the importance of inter- 
actions and disorder in the bath. Moreover, disordered 
spin baths warrant further attention because both inter- 
actions (often dipolar) and disorder are present in real 
spin bath systems like quantum dots in semiconducting 
hcteroj unctions and in Si:P. 
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In this paper, we re-examine the decoherence induced 
by disordered interacting spinbaths at finite tempera- 
tures. More precisely, we study the effect of an Ising 
bath with random spin-spin interactions on the coher- 
ence of a single spin 1/2. The random interactions are 
characterized by their variance A 2 where A is analogous 
to the cut off frequency for a bosonic bath as well as a 
mean value Jo which has no bosonic counterpart. Our 
choice of an Ising bath is primarily to facilitate an an- 
alytical study of the problem at finite temperatures in 
the thermodynamic limit. To ensure a quantum deco- 
herence of the central spin in our model, the central spin 
is coupled to the transverse spin components of the bath. 
To better comprehend the effect of the bath, we consider 
a model in which the time evolution of the central spin 
is exclusively governed by its coupling to the bath. For 
such a model, the problem is exactly soluble for a bath 
comprising independent spins/bosons. However, when 
intra-bath interactions are present, no exact solution can 
be obtained and one has to take recourse to approximate 
methods. In this paper, we only study the limit of a weak 
coupling of the central spin to the bath, where robust 
analytical methods are available to study the problem 
in an unbiased manner. We use the resolvent operator 
method, which takes us beyond the oft used Markovian 
master equation approaches to study the decoherence of 
the central spin induced by a perturbative coupling to 
the Ising spin bath. An interesting aspect of our work 
is that in the absence of any dynamics intrinsic to the 
central spin, for weak coupling to the bath, the decoher- 
ence is primarily dictated by the local-field distribution 
or equivalently, the dynamic structure factor of the bath 
spins. For a bath described by an Ising spin chain, we 
obtain the Markovian decoherence time scale and the non 
Markovian corrections as a function of the temperature 
and the parameters of the bath. We also discuss the case 
of Ising spins on various lattices in the high temperature 
regime and the case of the infinite-ranged Sherrington- 
Kirkpatrick spin glass model. 

The paper is organized as follows: we present the 
model and derive a general expression for the decoher- 
ence of the central spin in Sec. ITT1 followed by a discussion 
of the weak coupling regime in Sec. lIIII We then present 
our results for different disordered Ising spin baths in 
Sec. El 

II. MODEL 

We present the model used to study the decoherence of 
the central spin er c , weakly coupled to a bond-disordered 
bath of N Ising spins <Xi in the thermodynamic limit 
N — > oo. The total Hamiltonian describing the combined 
system of the central spin and the spin bath is given by 

H = H B +a x c V (1) 

= £^<^ ( 2 ) 

(ij) ' 



where, a% is the x-component Pauli operator of the cen- 
tral spin and of and of denote the Pauli operators of 
the bath spins and J^ are the interaction strengths be- 
tween the bath spins. Depending on the details of the 
model studied, (ij) could represent interactions between 
nearest neighbour spins or interactions of infinite range. 
In contrast to the models studied in Ref. where the 
bath hamiltonian Hb and the bath operator V that cou- 
ples to the central spin commute, here our choice of V is 
such that [Hb,V] ^ 0. The central spin and spin bath 
coupling is characterized by the parameters A^. Since, 
we are interested in the influence of disorder as well as 
the tendency of the system to order, we consider random 
interaction energies Jy which are quenched random vari- 
ables drawn from a distribution p(J) with mean Jo and 
variance A 2 . Though the central spin does not have any 
intrinsic dynamics, its coupling to the bath generates a 
non trivial dynamical behaviour. We note that Hb is 
the usual Ising spin glass Hamiltonian which has been 
well studied in the pasliii Depending on the distribu- 
tion of the spin-spin interactions and the dimensionality, 
this model can exhibit ferromagnetic, antiferromagnetic 
or even spin glass order in some temperature range. Since 
these phenomena have ramificiations for the collective be- 
haviour of the bath, it is reasonable to expect the result- 
ing decoherence to depend crucially on the underlying 
order in the bath. 

It is important to note that the formalism developed in 
this section and Sec. Illll is a priori applicable to any bath 
hamiltonian Hb (bosonic baths, Heisenberg spin baths, 
baths with both spins and bosons etc). For a Hamiltonian 
of the form (J2J, since a* is a constant of motion, it is 
convenient to directly study the time evolution of the 
reduced density matrix of the central spin 

p(t) = Tr B {e' lHt VLe lHt ) (3) 

where f2 is the initial density matrix of the composite 
system consisting of the central spin and the bath and 
Trs denotes the partial trace over the bath degrees of 
freedom. We use the units % = fcs = 1 in this pa- 
per. Denoting the eigenstates of a* by | <— ) and | — >), 
we see that due to the property of the density matrix 
Tr[p(i)] = 1 and the stationarity of Tr[p(t)a x ], the diag- 
onal elements of the density matrix given by the popu- 
lations (<— \p(t) \ «— ) and (— ► \p(t) \ — >) remain constant. 
Only the coherences (— ► \p(t)\ «— ) = («— \p(t) \ — ►)* rep- 
resenting the off-diagonal elements of the density matrix 
p(t) change with time. Often these coherences vanish at 
long times resulting in a decoherence of the central spin 
i.e., the asymptotic state of the central spin is a mixture 
of the states | <— ) and | — ») irrespective of the nature of 
the initial state f2. 

To simplify the calculation, we assume that at time t = 
0, the central spin and the bath are disentangled resulting 
in a factorizable initial density matrix: Q — p(0) ® pb- 
We further suppose that initially the central spin is in 
a pure state \ip) = a\ <— ) + (3\ — ►) (p(0) = \ip)(tp\) an d 
that the bath is at thermal equilibrium with temperature 
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T=l/0: 



Pb = 



,-0H B 

Z 



(4) 



where Z = Tr exp(— (3Hb) is the bath partition function. 
With these initial conditions we obtain the time evolved 
reduced density matrix 



P(t) 



-><< 



>)(->] + M(t)a*/3\->)(* 
I 



•| + M(t)VIHHI 



(•5) 



where the factor 



M(t) = Tr 



' e -i(fl-B+v)« ps e i(H B -v)t 



(6) 



is a measure of the decoherence induced by the bath at 
time t. Here Tr denotes the usual trace as Hb and V are 
operators in the bath Hilbert space. Under a rotation 
of 7r around the z-axis, while Hb and pb remain un- 
changed, V — > — V. Consequently the coherence M(i) is 
a real number. As mentioned in the introduction, though 
the coherence M is easy to evaluate for baths consisting 
of non- interacting spins/bosons (see Appendix |A"|) . it is 
rather difficult to estimate for an interacting spin bath for 
arbitrary values of the coupling to the central spin. We 
remark that M(t) is related to Loschmidt echoes of the 
bath which characterize its sensitivity to perturbations 
in its equations of motioniii£*i&. 



III. WEAK COUPLING REGIME 

In this section, we present a perturbative formalism to 
calculate M(t) valid for weak coupling to the bath i.e., 
the energy scale of the operator V is much smaller than 
all the energy scales of the bath. We use the resolvent 
operator method which goes beyond the Born-Markov 
approach and, though perturbative, is intrinsically ca- 
pable of handling non-Markovian time evolutions. We 
obtain a tractable expression for the decoherence in the 
weak coupling regime. 



A. Resolvent operator method 

To determine the complete time evolution of the co- 
herence M in the weak coupling limit V — ► 0, we first 
express it in terms of a self energy £iLi&. To obtain 
the self energy, it is convenient to work with the Laplace 
transform of M(t) 



M(z) 



dte lzt M(t) 



(7) 



where z is a complex variable with Imz > 0. As shown 
in Appendix IbI this Laplace transform can be written as 



where the self-energy £ is given by 

E(«) = Tr^CvPB +C v Q{z-QC Q^QCvPb] (9) 

In the above expression, the superoperators Q, Lb, C-v 
and C are defined by their actions on any operator A 
acting on the bath Hilbert space: QA = A — Tr(A)ps, 
C B A = [Hb, A], C V A = VA + AV and C = C B +C V . 
As is a density matrix, Q is a projection operator, 
i.e. Q 2 = Q. Note that Cy is not a Liouville operator 
whereas Cb is the Liouville operator corresponding to 
the bath Hamiltonian Hb- We reiterate that the above 
derivation for iJHJ and (0 is independent of the specific 
Hamiltonian Hb and coupling operator V considered in 
this paper. The self-energy £ can now be expanded per- 
turbatively in the interaction operator V. The second 
order result is given by the expression © with C re- 
placed by the bath Liouvillian Cb- The first-order term 
Tr(ZVps) = 2Tr(Ups) vanishes for the Ising spin bath 
Hamiltonian Hb and the interaction operator V defined 
by Therefore, the first non-zero contribution to the 
self-energy is given by the second-order term E2 which 
can be rewritten in terms of the time-dependent sym- 
metrised correlation function of V (see Appendix Q : 



E 2 (z) 



-2i 



dte lzt [(V(t)V) + (VV(t))} (10) 



M(z) 



£(*) 



(8) 



where (...) refers to the thermal average over bath spin 
configurations. Neglecting higher order contributions to 
E in |(HJ| is equivalent to the Born approximation 2 . We 
will see in the following that this approximation can de- 
scribe the decoherence at all timescales whereas a direct 
expansion of the coherence M gives only the short-time 
evolution. We remark that in the case of Heisenberg 
spins, an underlying magnetic order could result in a first 
order contribution to the self energy which would then 
lead to an oscillatory behaviour of M(t). 

The advantage of the resolvent operator formalism is 
that we can use the analyticity properties of the self en- 
ergy E to obtain a tractable expression for the coherence 
M. As shown in Appendix^! since Hb and V are Hermi- 
tian operators, the spectrum of the operator £ is real and 
hence the self-energy E is analytic in the upper (lower) 
half plane. Furthermore, since the spectrum of C in the 
thermodynamic limit is expected to be a continuum for 
the models considered in this paper, the self-energy E 
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manifests a branch cut on the real axis. The coherence 
M can thus be written in terms of the real functions A 
and r defined by 

A(E)-iT(E) = lim E(£ + ot) (11) 

where E is real. Since M(t) is real, the functions A and T 
satisfy A(-E) = —A(E) andT(-£0 = T(E). Performing 
the inverse Laplace transform of © and taking the limit 
?/^0 + , we obtain 

i f e~ itE 

e(t)M(t) = — / dE -— t=t (12) 

w w 2tt J E- A(E) + iT(E) v ; 

where Q(t) is the Heaviside step function. Moreover, 
integrating T,(z)/(z — E) along an appropriate contour 
in the upper half plane, one obtains the Kramers-Kronig 
like dispersion relation 

A(E) - iT(E) = - l -V J [HE') iT(E')\ (13) 

where V denotes the Cauchy principal value. This shows 
that Q(t)M(t) is exclusively determined by T (or A). 



B. Non-Markovian evolution 

In this section, we analyze (I12fl to determine M(t) at 
any time t in the weak coupling regime. We first note that 
differentiating (|12J1 with respect to the time t yields the 
conditions M(0) = 1 and d t M(Q) = 0. Differentiating 
(|12fl further and taking the limit V — ► 0, we find, to 
second order in V, 

Q(t)d 2 M ~ -J- [ dEe- iEt [T 2 (E)+iA 2 {E)} 
2ir J 

- J dEe- lEt T 2 (E) (14) 

where T 2 and A2 denote the second-order terms of the 
functions T and A. Note that to obtain the second equal- 
ity we have used the relation l|13|l and the Fourier trans- 
form of 28 (t) — 1. A solution to the above differential 
equation, for t > 0, is 

«(«),i-?/ ffi ^r, (S , (15) 

where we have taken into account the symmetry of the 
function T 2 . Though this Fermi golden rule like approx- 
imation yields the correct short time behaviour, it leads 
to the false result M(t) ~ —T 2 (0)t for t — > 00 and hence 
is invalid for an arbitrary time t. 

To obtain the long-time decoherence, we evaluate the 
integral (|12fl using the analytic continuation of S from 
the upper half plane to the lower half plane (second Ric- 
mann sheet). If the functions T and A are analytic in the 



vicinity of E — 0, the analytic continuation of E in the 
second Riemann sheet is given by 

t(z) = -iT(0)+A'(0)z + O(\z\ 2 ) (16) 

for small z where A' denotes the derivative of the func- 
tion A with respect to E. Note that the symmetry prop- 
erties of the functions A and T ensure A(0) = and 
r'(0) = 0. The coherence ()12f) is principally determined 
by the singularities of [z — E(z)] -1 , one of which is a 
pole at zq = —iT 2 (0) + 0(V 4 ) close to the real axis. 
Every other singularity lies beyond a finite distance e _1 
from the real axis determined by the temperature and 
the bath parameters. Consequently, for times f>e, the 
residue exp(— izot)[l — 9 2 E(zo)] _1 dominates and thus in 
the weak coupling limit V — > 0, 

lnM(t) ~ -r 2 (0)t + A£(0) . (17) 

Combining the approximations for short times and long 
times given by l|15|l and (|17f) respectively, we see that in 
the weak coupling regime the decoherence at any time t 
is described by 

bM^-f/^fi^r^) . (i 8 , 

This equation shows that M(t) in the weak coupling 
regime is determined by the entire function T 2 . For 
asymptotic times, we see from (I18f) that the coher- 
ence of the central spin is essentially given by M(t) ~ 
exp[— r2(0)t] which is simply the solution of the Marko- 
vian master equation obtained within the Born-Markov 
approximation 2 . However, for the short and interme- 
diate time evolution of M, the full energy dependence 
of T 2 comes into play which can then lead to a non- 
exponential decay i.e., non-Markovian behaviour of the 
coherence. Depending on the temperature and bath pa- 
rameters, the asymptotic Markovian regime could even 
disappear, provided that A 2 (0) goes to infinity. We re- 
mark that Ea. l|18|l is valid to all orders in V for a bath 
of independent bosons. 

IV. DISORDERED ISING SPIN BATHS 

In this section, we use the formalism of the previous 
section to determine the decoherence induced by vari- 
ous disordered Ising spin baths. Though the spinbath 
models of real systems are expected to be more compli- 
cated than the Ising Hamiltonians considered here, we 
nonetheless study these systems in various dimensions 
to understand the effect of these simpler systems on the 
decoherence. More precisely we study the effect of the 
disordered Ising spin chain on the coherence M of the 
central spin as a function of the temperature and the 
bath parameters Jo and A. We also make predictions 
for the infinite-ranged Ising model also known as the 
Sherrington-Kirkpatrick model and for Ising models in 
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higher dimensions in the high temperature regime. Be- 
fore we embark on in-depth calculations of M, we show 
that for any Ising bath, the decoherence of the central 
spin in the weak coupling regime is intimately linked to 
the local-field distribution of the bath. 



A. Local-field distribution 

As shown in Sec. IIHI the decoherence in the weak 
coupling regime is determined by the time-dependent 
correlation function Re(V(t)V), through iJTSJ and (|TT)j) . 
For the case of the Ising bath Hamiltonian defined in 
(J2J), this correlation can be obtained from the local- 
field distribution of the bath. To see this, it is 
useful to work in the eigenbasis |{<7i}) of Hb ■ 
Hb\{<Ji}) = — ^u' f7 » cr j|{ (7 »})- As tne matrix ele- 
ment ({ci}|crf |{(Tj'}) is nonvanishing only for spin con- 
figurations {cr,} and {o[} where a[ = <7j(l — 25^), the 
time-dependent correlation of V is a sum of local corre- 
rations : (V(t)V) = £ fe A|<of We find 

Re{V(t)V) = \Y. a J^ Em J ' j<T,CTj cos [ 2t YlT Jki(7i 

(19) 

where Y]j denotes a sum over the spins <Tj interact- 
ing with the spin cr k . Note that the term J3i Jki^i 
in (|19|l is the effective local field acting on the spin 
at site k generated by the configuration {at}. Since 
Re(a x (t)a x ) = cos(2th) for an isolated spin <r in a field 
h parallel to the z-axis, we rewrite <|19ll as 

Re{V(t)V) = Y X l I dh p k{h) cos(2th) (20) 
k J 

where P k {h) = (^8 (h — Jki^i^j ^ can be interpreted 
as the distribution of the local field h at site k at tem- 
perature T. This interpretation is also valid for the bath 
thermodynamic quantities such as magnetisation or spe- 
cific hea1»i^. 

It is now rather straightforward to obtain the func- 
tion T2 which is the crucial ingredient to determine the 
coherence M in the weak coupling regime. To avoid com- 
putational complexity, in the rest of the paper, we con- 
sider equal couplings \ k — AA -1 / 2 . With this choice, 
the function reads 

r 2 (E) = 2tt\ 2 P(E/2) (21) 

where the local-field distribution P is the spatial average 

The thermodynamic limit of this expression is unambigu- 
ously defined and P is self-averaging^ i.e., in the limit 
A — > oo, P is given by bond disorder average of any 



distribution P k . For an Ising spin system, the local-field 
distribution P is temperature dependent and determines 
the thermodynamic quantities and the dynamic linear 
response of the systemi 9-21 . Eq. 121|) shows that the 
local-field distribution is also an important characteristic 
of an Ising system considered as a bath. 

We remark that for a bath of independent spins-, i.e. 
H B = — hio\ in @, the function T 2 is given by J2D 
with the field distribution P(h) = J2k( X k/ 'X) 2 S(h - h k ). 
In this case, since the spins are non-interacting, the distri- 
bution P is arbitrary and temperature independent and 
the ensuing decoherence of the central spin is tempera- 
ture independent. This feature of a temperature indepen- 
dent decoherence induced by a bath of independent spins 
seen in the weak coupling limit is also seen in the exact 
non perturbative result for the coherence M(t) obtained 
in Appendix [a1 

B. Ising spin chain 

In this section we consider a ID Ising bath described 
by the Hamiltonian 

N-1 

H B = -Y J Ji<°i + i ■ (23) 
»=i 

The spin at site i interacts with its nearest neighbors 
with interaction strengths Ji and Ji-\. The Jj are a 
quenched set of random bonds drawn independently from 
a distribution p with mean Jo and variance A 2 . 

1. Born self-energy 

To obtain the coherence M we need the time- 
dependent correlation of the coupling operator V — 
~^2i X i <T i ■ -As shown earlier, this correlation is given 
by Re(V(t)V) = E fe A fc Re K W^)- Here the time de- 
pendent spin-spin correlations can be written in terms of 
the static correlation function as 

Re(a x (t)a%) = cos(2t J fe _i) cos(2i J k ) (24) 
-(o-fe+iOfc-i) sin(2<J /£ _i) sin(2iJ fe ) 

where (&k+i a k-i) = tanh(/3Jfc_i) tanh(/3Jfc) is related to 
the derivative of the partition function 

N-1 

Z = 2N II C0S HPJi) (25) 

i=l 

with respect to the interaction strengths J k and Jk-i- 
The choice X k = AA" 1 / 2 yields a variance of Re(V(t)V) 
of order A -1 and a mean of order A°. Consequently, 
Re(V(t)V) and hence £2 are self-averaging in the ther- 
modynamic limit. The second-order self-energy is thus 
given by the average of l|10|) over bond disorder. We 
obtain 
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T 2 (E) = 2ir\ 2 / dJ Pe {J) Pe {J + E/2) [1 - tanh(/3J) tanh(/3 J + (3E/2)} 



(26) 





= 


T = 0.02 A 


J t 


= 


T - A 


J t 


= 


T= 10A 


J t 


= 2A 


T = 0.02 A 




FIG. 1: The real part of the dimensionless second-order self- 
energy F2A/A 2 as a function of E/A for a Gaussian bond 
distribution with mean values Jo = 0, 2 A and temperatures 
T = 0.02A, A, 10A. 



where p e (J) = [p{J) J)]/% is the symmetrized bond 
distribution. Note that the integral on the right side is 
the local-field distribution (|2*2)> for the spin chain 20 ! 22 . 
Equation 126(1 is valid for all bond distributions p. More- 
over, for distributions p symmetric around their mean 
values, as T2 depends only on the symmetrized distri- 
bution p e the self-energy £2 is the same for opposite 
means ±Jo- This shows that, in this case, though the 
decoherence is influenced by the interactions in the bath, 
it is insensitive to the underlying ferromagnetic or an- 
tiferromagnetic nature of the interactions. This feature 
can be understood as follows. The spin chain Hamil- 
tonian (|23[1 and the local field operator Jfe-ierjjLj + 
Jk&k+i are invariant under the transformation [J,-, of] — > 
[— Ji, (— l) l_fe of]. Consequently, the local-field distribu- 
tion P and hence the function T2 are invariant under the 
transformation Jo — > — Jo. To complete the determina- 
tion of the second-order self-energy, its real part A2 can 
be obtained from the function T2 using i|13|) . In Fig. ^ 
we plot T2 for different values of Jo and the temperature 
T. 



2. Markovian evolution 



As shown in Sec. lIIII in the weak coupling limit A — > 0, 
the evolution of the coherence M is essentially Markovian 
: M(t) ~ exp(— jt). For the disordered Ising chain, the 



decoherence rate 7 is given, to lowest order in A, by 

7 = r 2 (0) = 2ttA 2 J d.Jp e (J) 2 [l - tanh(/?J) 2 ] . 

(27) 

The rate 7 increases monotonically from 7 = at T = 
and saturates to the value 



7o 



27rA 2 y d.Jp e (J) 2 



(28) 



as T — > 00. At low temperatures, 7 ~ 47rA 2 p(0) 2 T for 
disorder distributions with p(0) 7^ 0. The vanishing of 7 
at zero temperature can be understood as follows 2 ^. At 
T = 0, the local field distribution is completely dictated 
by the spin structure of the ground state |{ci}) of the 
bath Hamiltonian (|23|l . Since JiOiOi+x > for any pair 
of neighboring spins in the ground state, this inevitably 
leads to nonvanishing local fields Ji-i<7i_i + Ji<Ji+\ = 
{\Ji-i\ + l-ZiDci- This implies P(h = 0) = and hence 
7 = 2ttA 2 p(0) = 0. 

We now study the influence of the bond distribution 
on the rate 7. In the weak disorder regime A <C | Jrj|, the 
distribution p e ( J) 2 practically vanishes for J 7^ ±Jo thus 
the rate 7 is given by 



7 ~ ttA 2 / dJp(J) 2 [1 - tanh(/3J ) 5 



(29) 



In this regime, the temperature variations of 7 are inde- 
pendent of the form of the bond distribution p and are 
exclusively determined by the mean value Jo- In Fig. [21 
we plot the temperature dependence of 7 for various val- 
ues of the mean interaction strength Jo in the case of a 
Gaussian bond distribution 



; -(./-./o) 2 /2A 2 

V2^A 



(30) 



The linear regime at low temperature given by 7 ~ 
2(A/A) 2 exp(-J 2 /A 2 )T is visible only for |J | < 2A. 
For higher interaction strengths Jo , 7 remains practically 
zero in the low-temperature regime. The maximal rate 
obtained as T — * 00 is given by (|2*8)l : 700 = ^/ttX 2 [1 + 
exp(- J 2 /A 2 )]/2A. For | J | > 2A, the high-temperature 
rate 700 is essentially independent of Jo and the tem- 
perature dependence of 7 is well described by the weak 
disorder approximation 7 ~ (y / 7rA 2 /2A)[l — tanh(/3J ) 2 ]. 
Note that the agreement with this approximation im- 
proves with increasing Jo. 

We now consider a uniform distribution for the intra- 
bath interaction strength: 



p{J) = (2V3AY 



for |J-J |<a/3A 
otherwise . (31) 



7 



In this case, the integral l|27[l can be evaluated exactly 
I 




FIG. 2: Dimensionless rate 7 A/ A 2 as a function of the tem- 
perature T/A for a Gaussian bond distribution with mean 
Jo = 0, A and 2A. The weak disorder approximation 129^ is 
shown for Jo = 2 A. 

For I J 1 < V3A, 7 oc T at low temperatures with a slope 
7rA 2 /3A 2 independent of Jo whereas, for | J | > \/3A the 
low temperature behavior is not linear. For | Jq\ > 4A, 
the agreement between the exact result and the weak dis- 
order approximation 7 ~ (7rA 2 /2\/3A)[l — tanh(/3J ) 2 ] 
is excellent. Moreover, in this regime, the Gaussian 
and uniform bond distributions cannot be distinguished 
~ 1.02). An interesting feature of our results is 
that a non-zero Jo favours the coherence of the central 
spin via a robust short range ordering of the bath spins. 

Comparing the above results with those obtained for a 
bath comprising free spins, we see that the timescale of 
the decoherence generated by the spin chain bath 7^ ~ 
A /A 2 (for weak coupling A <C A) is much longer than 
the decoherence time A -1 (for A <C \/~N) obtained in 
the free case (see Appendix lATl . As we will demonstrate 
later, 7^, 3> A -1 for any dimensionality of the spin bath 
lattice. This clearly illustrates the fact that interactions 
in the bath significantly slow down the decoherence of 
the central spin. 

Another interesting comparison is to an Ohmic bo- 
son bath. Contrary to the Ising bath, the decoherence 
rate jbos of an Ohmic boson bath is proportional to T 



for any temperature T : 



(32) 



I 

in the whole temperature range and thus does not sat- 
urate at high temperatures. This forces the question as 
to whether interactions between the bosons also lead to 
a saturation of the rate 7(, os . This warrants further work 
which is beyond the scope of the present paper. If the 
central spin is coupled both to a spin bath and a boson 
bath, the resulting coherence M(t) is given by the prod- 
uct of © and a similar factor, with Hb replaced by the 
boson bath Hamiltonian Hb and the interaction opera- 
tor V by an analogous boson operator V&. Then in the 
Markovian regime at weak coupling, the total decoher- 
ence rate is simply the sum jbos + 7- Consequently, at 
high temperatures the Ohmic boson bath dominates but 
for temperature T<A the Ising bath has to be taken 
into account. The question of their relative dominance 
depends on the various bath coupling constants and can 
vary from system to system. 



3. N on- Markovian regime 

Here, we discuss the non-Markovian aspects of the de- 
coherence of the central spin essentially seen at low tem- 
peratures and at short and intermediate time scales. We 
first consider the case of a bond distribution with zero 
mean. The time evolution of the coherence M(t) within 
the Born approximation, shown in Fig. |3J is obtained by 
a numerical evaluation of T2(E), A2(E) and M(t) using 
the expressions (|26|) . <|13fl and 112|) for a Gaussian bond 
distribution l|3U|l with Jo = 0. Note that the agreement 
with the weak coupling approximation i|18|) is remarkably 
good even for reasonably large values of A i.e., of the or- 
der of 0.1A. Though, in Fig. [3 we illustrated the equiv- 
alence of the Born and weak coupling approximations for 
M(t) in the weak coupling regime, we nonetheless expect, 
based on the analyticity arguments presented in Sec. II I II 
the weak coupling approximation i|18[l to be an exact de- 
scription of the full coherence M for low enough A. 

We now discuss the influence of temperature on the 
behavior of the coherence M(t). In Fig. [31 we see that 
at high temperatures, M remains practically constant 
for times t< A -1 and then decays exponentially as de- 
scribed by l|17|l . For temperatures T<A, the Marko- 
vian regime is preceded by a novel intermediate time 
regime A _1 <i</3. The difference between the high- 



3 tanh 



tanh 



|Jo|-\/3A 
T 

Jo I - V3A 
T 



for I J 1 < V3A 
for |J |>V3A 



temperature and the low-temperature decoherence can 
be traced back to the temperature dependence of the 
function T 2 . At high enough temperatures, since T 2 (E) 
is essentially a peak of width A, one crosses over from the 
short time regime to the Markovian regime for t ~ A -1 . 
As the temperature is lowered, ^(0) steadily decreases 
but the curvature r 2 (0) remains negative. However, be- 
low a certain temperature r' 2 '(0) becomes positive, see 
Fig. ^ an d the function T 2 can be effectively character- 
ized by two energy scales: TandA. At low temperatures, 
T 2 remains practically constant for \E\<T, increases 
linearly for larger energies with a slope (A/A) 2 and fi- 
nally vanishes for \E\>A. These three energy regimes 
result in three different time regimes for the coherence 
Q18JI. The short-time behaviour (t< A" 1 ) is determined 
by the long-energy tails of T 2 . For intermediate times 
^t<(3, the linear regime of T 2 yields a power law 
decay of the coherence M(t) oc t~ e where the exponent 
is given by 



For long times (t > /?) the integral H18|l is dominated 
by the energies \E\<T for which T 2 (E) ~ 7 and hence 
the decoherence is essentially exponential with the rate 
7. However, one should be cautious about extending the 
above results to ultra-low temperatures because the con- 
tribution of the energies \E\ > T, given by A 2 (0), diverges 
in the limit T — > 0. This divergence is logarithmic with a 
prefactor e. At zero temperature, the Markovian regime 
disappears and the coherence vanishes in the limit t — > 00 
according to the power law t~ e . The low temperature 
behavior of M(t) obtained here is very similar to the de- 
coherence induced by a boson bath in the strict Ohmic 
case^ i.e. for an Ohmic spectral function with a cut-off 
frequency w c far larger than T and 1/t. In this case, the 
coherence of the spin coupled to the boson bath is given 
by In M(t) = -K \u[(uj c /Ttt) smh(TrTt)} where K is the 
spin-bath coupling strength. We recover a Markovian be- 
havior, \nM(t) ~ K \n(2TrT/u; c ) - KnTt, for times t > (3 
and a power law, M(t) = (uj c ty K , at zero temperature. 

We now present a more detailed comparison of our 
numerical results with the weak coupling approximation 
(|18ll . shown in Fig. [3] For given A and T, the Born ap- 
proximation deviates from the expression (|18|l as time 
increases. This difference, which is more pronounced at 
low temperatures, can be easily quantified in terms of 
higher order corrections stemming from the expansion of 
the analytically continued self-energy discussed in Sec. 
IIII Bl At long times, the first correction to lnM(t) arises 
at 0(A 4 ) and takes the form 

InM(t)-[A4(Q)-r 2 (0)i] ~ -r 2 (0)r' 2 '(0) + ±[A 2 (0)] 2 

-r 2 (0)A 2 (0)t . (34) 

In the high temperature regime, since A 2 (0) is positive, 
the weak coupling approximation given by the second- 
order terms in l|34|) slightly overestimates InM. As the 




tA 



FIG. 3: (A/A) 2 In M as a function of time tA within the 
Born approximation for a Gaussian bond distribution with 
zero mean for coupling strengths A = 0.1A, 0.2A and tem- 
peratures (a) T = 10A and (b) T = 0.02A. The full line is the 
weak coupling approximation and the dotted line corresponds 
to a free spin bath. 



temperature is lowered, the evolution of the function T 2 
discussed earlier (see Fig. ^) results in a sign change 
of A' 2 (0). Thus for low enough temperatures, A' 2 (0) is 
negative and the weak coupling approximation under- 
estimates In M. This analysis is indeed very consistent 
with our results shown in Fig. [3] This discussion clearly 
highlights the efficiency of our approach based on the 
analyticity properties of the self-energy, obtained within 
any approximation, to describe the long-time decoher- 
ence. We remark that other fourth order corrections to 
]uM(t) exist beyond the Born approximation. 

We now consider the case Jo 7^ 0. For weak disorder 
I Jo I ^> A, the function T 2 at low temperature consists 
of two Gaussian peaks of width 2^/2A centered around 
±4 Jo. The low-tcmpcrature self-energy in this regime 
is thus very different from that for Jo = 0, see Fig. ^ 
Nonetheless, the qualitative behavior of the coherence 
M at intermediate and long times, determined by the 
low-energy T 2 , is very similar to the one discussed above 
for Jq = 0. For short times t< A -1 , contrary to the 
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lnM(t) oc — t 2 behavior seen for Jo = 0, here the coher- 
ence M(t) oscillates at the frequency 2J /tt. At zero tem- 
perature, the long-time decoherence decays as a power 
law, M(t) oc t~ v where the exponent r\ = eexp(— Jq/A 2 ). 
As seen earlier for the Markovian decay rate 7, the low- 
temperature behavior seen here also signals a slowing 
down of the decoherence by a nonzero Jo- 

C. Higher-dimensional spin lattices 

Here we consider other geometries for an Ising 
spin bath described by the Hamiltonian Hb = 
— Yluj) Jij cr i' T j where the spins occupy the sites i of a 
regular lattice of arbitrary dimensionality and are cou- 
pled by nearest-neighbor interactions. The bonds Jy are 



drawn independently from the Gaussian distribution p 
with mean J and variance A 2 given by (|3UI) . An im- 
portant difference between the spin chain model and the 
generic Ising model on higher dimensional lattices is the 
presence of frustration arising from geometric constraints 
and/or randomness. Though frustration can give rise to 
novel ground states and related dynamical behaviour, it 
renders any analytic study of these models very difficult. 
In this section, we focus on the high temperature regime 
where one can use a controlled analytical method like the 
high temperature series expansion to study the effect of 
higher dimensions on the decoherence. 

To obtain the coherence M at high temperatures, we 
expand the second-order self-energy IjlUII in terms of the 
inverse temperature (3. To do so, we first rewrite the 
time-dependent correlation (|19|l as 



Re(V(t)V) = — x t il^ 1 + £r « 7 j' s w) Re II j [ cos ( 2 ^fc) + sin(2tj ife )] (35) 

k (ij) 

I 



where J}^ ' denotes a product over the nearest neigh- 
bors of site k, Z' = J2{a,} n<ij)(! + and = 
tanh(/?Jy). Since Ky — » as (3 — > 0, we consider the 
above expression as a power series in Km . Multiplying out 
the two products in Ea. l|33|l generates a series of prod- 
ucts of nearest-neigbour spin pairs : {uiUjUktri...). Since 
Ea. (|35p involves sums over configurations {ci}, each of 
these spin pair products contributes only if it simplifies 
to 1. This implies, for instance, that the n th order terms 
in the expansion of Z' in Kij correspond to closed loops 
comprising n bonds on the lattice. An immediate conse- 
quence is that Z 1 = 1 up to third order. However, in the 
expression (|35|l there exists another type of relevant spin 
pair products which involve repeated bonds. We discuss 
these terms in the following. 

As T — > 00, for equal couplings to the central spin 
\k = XN^ 1 / 2 , the time-dependent correlation l|35|l be- 
comes Re(V(t)V) = X 2 N~ 1 J2k Tli^ cos(2i Jj fc ). Note 
that this correlation remains unchanged under J^ — > 
— Jik- This invariance stems from the equiprobability 
of every spin configuration at infinite temperature. This 
infinite-temperature correlation is self-averaging in the 
thermodynamic limit and hence the corresponding sec- 
ond order self-energy is given by the average of (|f Oj l over 
the bond distribution. This leads to 

/oc 
dt e lEt e~ 2A2t2s cos(2 J t) s (36) 
-00 

where s is the coordination number of the lattice. As 
remarked earlier for the ID case, here also I^ is the same 
for ±Jo for bond distributions symmetric around their 
mean values. Eq. I|36|l shows that in the infinite tempera- 
ture limit the only characteristic of the bath lattice which 



intervenes in the decoherence is s. In particular, there is 
no explicit dependence on the dimensionality. However, 
we shall show later that the detailed geometric charac- 
teristics of the lattice manifest themselves in the higher 
order corrections. Though static properties of the spin 
system are independent of Jo and A in the infinite tem- 
perature limit, these parameters strongly influence the 
correlation ~Re(V(t)V) and hence the decoherence. Us- 
ing H36fl . for the case Jo = 0, we obtain from (|f 8|) 

lnM(f) = -2V^A 2 -rterf(t/r) + 2A 2 t 2 (l - e^ 2 ^) 

(. 3? ) 

where erf is the error function and the characteristic time 
r is defined by 

t = A- 1 (2s)- 1 / 2 . (38) 

For times f>r, we recover the Markovian regime where 
lnM(t) ~ ~2^/n\ 2 Tt + 2A 2 r 2 . For t < r, one finds the 
usual short time evolution : lnM(t) ~ — 2X 2 t 2 . The same 
time regimes exist for Jo 7^ 0. For short times i « t we 
find, for an even s, 

A 2 ^ slsm[J t(2n~ S )] 2 s\X 2 t 2 
11 2S " 2J o h'o n\{s - n)\{2n - s) 2 2 s -![(s/2)!] 2 

If s is an odd number, the short-time decoherence is de- 
scribed by (|39|l without the quadratic term in t and with 
a sum running from to (s — l)/2. Clearly the effect 
of a nonzero J is to induce oscillations in the coherence 
M(t). The decoherence rate in the Markovian regime is 
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given by 



A 2 



2 s " 1 A V 2s 



i 

-y 



-Jg{2n-sf/2sA 2 



— ' nl(s — n)\ 



(40) 

For s = 2 this expression simplifies to the result i|28|) 
obtained for the spin chain. As anticipated, the deco- 
herence time j^f is of the order of A /A 2 and hence far 
longer than that for free spins which is of the order of 
A . The infinite-temperature rate manifests the influ- 
ence of the lattice geometry: for lattices with an odd s 



like the honeycomb lattice (s — 3), 7^ — > in the weak 
disorder limit A <C | Jq\ whereas for an even s, triangular 
and square lattices for example, 700 7^ in this limit. 

We now evaluate the leading order corrections in (3 
to the decoherence rate 7. To illustrate the significance 
of the lattice geometry we consider various bidimensional 
lattices. For a triangular lattice, the lowest order terms in 
the sum over the configurations in (|35[) take the form 
— sin(2< Jik)nij sin(2i Jjfc) where the site i is a neighbor of 
the site k and the site j is a neighbor of both sites k and 
i. Consequently, the rate up to first order reads 



J 



T 

Too 



A 2 /3 Jp 

~16 A - 



2 _ 2e- 4J '/ 3A2 + e~ J °/ 3A2 



-3Jq/ a 



0(/? 2 ) 



(41) 



where 7^ is given by (|40|l with s — 6. The explicit ap- 
pearence of Jo in this expansion is linked to the fact that 
the first nonzero correction for this lattice occurs at first 
order in /3. Moreover, since this correction is positive 
(negative) for baths with a majority of antiferromagnetic 
(ferromagnetic) bonds, the ensuing decoherence time is 
longer for a predominantly ferromagnetic bath. This 
clearly highlights the importance of frustration in de- 
termining decoherence. Thus at high temperatures, the 
decoherence induced by the triangular lattice is very dif- 
ferent than that by the linear chain : the convergence of 



7 T to 7^ is slower and depends on the sign of J . We also 
remark that though the cubic and triangular lattices have 
the same coordination s — 6 and hence the same infinite- 
temperature rate, the corrections are different since there 
is no three-bond closed loop on the cubic lattice. For the 
honeycomb lattice, the first corrections arise at second 
order. These corrections correspond to loops with two 
repeated bonds : —n ki sin(2tJ ik )n k jSn\(2tJjk). Here the 
contribution of closed loops becomes relevant at higher 
orders. The resulting decoherence rate is 



J 



IT X 2 



2 R 2 r 



24 A 



(3A 2 + 4J 2 ) e" J ° 2 / 6A2 + A 2 e -3./ 2 /2A 2 
I 



0(/3 4 



(42) 



where 7^ is given by (|4D|> with s = 3. Finally for the square lattice we obtain the decoherence rate 
I 



^1^! g A 2 - 2J 2 ) e - 2j2 / A2 + 3(A 2 + J 2 )e- j2 /2A 2 + ^A 2 + 



5Jn' 



o(p 4 ) 



(43) 



r 



where 7^ is given by l|40|l with s = 4. We mention 
that for the square lattice both closed loops and loops 
with repeated bonds contribute to the lowest order cor- 
rection. Interestingly, the convergence of the rates j H 
and 7 s to their respective infinite-temperature limits is 
reminiscent of that of the chain. Moreover, as seen for 
the spin chain, our results for ^ H and 7 s are independent 
of the sign of Jo. This feature can be attributed to the 



bipartite nature of these lattices. As explained for the 
chain, the respective local field operators and Hamilto- 
nians are invariant under a sign change of the bonds Jij 
coupled with a suitable unitary transformation for the 
spins. Since this argument is valid in the entire para- 
magnetic phase, the higher order corrections to j H and 
7 s are expected to be independent of the sign of Jo in 
this phase. On the other hand, due to its non-bipartite 
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nature no analoguous transformation exists for the trian- 
gular lattice. It would be interesting to study the evolu- 
tion of these rates behave as one lowers temperature and 
enters a non-paramagnetic phase. 



D. Infinite-ranged Ising bath 

We now consider a bath described by the mean field 
Sherrington-Kirkpatrick model where all the Ising spins 
interact with each other. We restrict ourselves to the well 
studied case of spin-spin interaction strengths distributed 
around a zero mean with a variance A 2 . This model ex- 
hibits a spin glass phase transition at a finite temperature 
T sg = A below which the spins freeze 2 ^. The spin glass 
phase is characterized by the large number of metastable 
states present which then lead to anomalous dynamic 
behaviour. An interesting question is whether the deco- 
herence manifests novel features as goes from the high 
temperature paramagnetic phase to the low temperature 
spin glass phase. 

As shown in Sec. II VI the coherence M(t) of the central 
spin in the weak coupling regime can be obtained through 
a knowledge of the local- field distribution P(h). Though 
the infinite-range model has been extensively studied in 
the past using the replica approach and numerical simu- 
lations, it is not easy to obtain the local-field distribution 
for all fields and temperatures. Here, we do not delve into 
the problem of calculating P(h) but use existing results^ 
to make predictions for the coherence of the central spin 
coupled to this mean field bath. In the high temperature 
paramagnetic phase, the local-field distribution is 



P(h) = 



1 



2AV2^ 



e -(h-f3A 2 ) 2 /2A 2 + e -{h+l3A 2 ) 2 / 



2A- 



In the spin glass phase, calculations based on the replica 
formalism yield the following result for P(0) provided the 
temperature T ~ T sg : 



P(0) 



1 



1-1- 



T 



-5"" 



T 

t7 



These analytic results are sufficient to determine the de- 
coherence in the Markovian regime : the decoherence rate 
is given by 7 = v^A 2 exp[-i(T ss /T) 2 ]/A for T > T sg 

and 7 = v^/iA 2 [1 - (1 - T/T sg ) - ± (1 - T/T sg ) 4 } /A for 
temperatures in the vicinity of the spin glass transition 
temperature i.e., T ~ T sg . Clearly, in the Markovian 
regime, one does not see any sign of the spin glass tran- 
sition. Note that 7 saturates at infinite temperature to 
a value comparable to that obtained earlier for finite- 
ranged lattice models. This implies that even in the case 
of a maximally frustrated bath, the central spin decoheres 
at timescales longer than those for free spins. 

For lower temperatures, only numerical solutions exist 
for the local-field distribution22*2i. These results suggest 
a continuous variation of P(h) with temperature. The 



only significant signature of the transition is a flattening 
of P at T = T sg which has no manifest effect on the 
decoherence. Moreover, numerical extrapolations of P(0) 
to low temperatures indicate a rate 7 ~ (A/A) 2 T which 
is very similar to the linear T behaviour seen in the Ising 
spin chain system at low enough temperatures. Again, 
since 7 = at zero temperature, the decoherence is no 
longer Markovian. The form of M(t) is then dictated 
by the low-energy behavior of T2 through (|18|) . Based 
on the numerical inference P(h) oc h for T — 0, the 
central spin is expected to decohere as a power law at 
T = 20 - 21 . To conclude, we see that both the thermal 
transition and the spin glass order of the Sherrington- 
Kirkpatrick model do not have any palpable effect on 
the asymptotic decoherence in the weak coupling regime. 



V. CONCLUSIONS 

In this paper, we have studied the decoherence induced 
by a Ising spin bath with random intra-bath interactions. 
The resolvent operator formalism was used to determine 
the coherence M of the central spin for weak coupling to 
the bath. We then obtained detailed analytical results for 
the disordered Ising spin chain bath for arbitrary temper- 
ature. The decoherence was found to be independent of 
the sign of the mean intra-bath interaction strength Jo 
for symmetric bond distributions. Three regimes were 
identified in the time evolution of M(t) : a short time 
Gaussian decay an intermediate time power law behavior 
and the usual asymptotic Markovian regime. The relative 
sizes of these regimes are fixed by temperature. At zero 
temperature, the Markovian regime was found to vanish 
and the decoherence is essentially described by a power 
law decay. We also studied the decoherence induced by 
an infinite-ranged Ising spin glass bath and Ising models 
on lattices in dimensions greater than one. For all these 
baths, the Markovian rate was found to saturate to a fi- 
nite value at infinite temperature, which is much smaller 
than the corresponding rate for a free spin bath. Our 
results clearly indicate that intra-bath interactions sig- 
nificantly increase the timescales over which the central 
spin decoheres. 

For the infinite-ranged Ising spin glass, our conclusions 
based on existing results suggest that the thermal spin 
glass transition has no visible effect on the decoherence. 
Plausibly this is an artefact of the infinite-ranged inter- 
actions and/or the Ising nature of the spins. This raises 
the general issue of the influence of thermal and quan- 
tum phase transitions and the resulting orders in finite- 
ranged spin baths on the decoherence of the central spin. 
In most realistic cases, the spin environment consists of 
Heisenberg spins. In this case, one expects the dynamics 
of the bath to be richer and this may have interesting 
consequences for the decoherence. This however, is be- 
yond the scope of the analytic work presented in this 
paper. An interesting question is the effect of a strong 
coupling between the central spin and the bath. For a 
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bath of independent spins/bosons the results obtained in 
the weak coupling regime are qualitatively valid even for 
strong coupling. However, in the presence of interactions, 
this is not necessarily the case and one can expect novel 
dynamical behaviors. A natural extension of our work 
would be to include the intrinsic dynamics of the central 
spin and study the relaxation induced by the spin bath. 
These and other questions are left for future work. 
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APPENDIX A: INDEPENDENT SPIN BATH 

Here we derive an exact expression for the coherence 
M(t) given by JSJl in the case of a bath comprising non- 
interacting spins described by the Hamiltonian Hb = 
— J2i^ l i (J i- Since the spins are independent, the trace 
© can be factorized as 



JY 



M(t) = n 



Tr [< 



,0h k al -it(h k a k -\ k al) it(h k a k +\ k a k )] 



k=l 



2cosh(/3/i/ £ ) 



(Al) 

To evaluate the factors in this expression we require 
the diagonal elements of 2 x 2 matrices of the form 
cxp[— ibA(a)) exp[ibA(a)] where A(a) — acr z + a x . Di- 
agonalizing A(a), we find that these diagonal elements 
are equal and given by 



(a\e~ tbA{a) e lbA(a) \<T) = 



1 



1 + a 2 



cos 



1 



(A2) 

The resulting coherence M(t) is independent of the tem- 
perature and can be written as 



N 



lnM(i) = J2 ln 



k=l 



1 - 



2A 2 



(A3) 

For hk = 0, the above expression leads to the coherence, 
M(t) = Ylk cos(2Afet) which then culminates in a Gaus- 
sian decay M(t) — exp(— 2i 2 ^ fe Af) for weak coupling 
to the bath. For nonzero local fields hk the coherence in 
the weak coupling regime reads: 



lnM(<) = -2A 2 / dh P(h) 



sm(ht)" 



h 2 



(A4) 



where the coupling strength A and the field distribution P 
are defined by A 2 = £ fe A 2 and P(h) = E^/A) 2 ^ - 
hk). This result is the same as that obtained earlier in 
(1 1 HI) for the specific bath considered here. 



APPENDIX B: DERIVATION OF THE 
SELF-ENERGY 

In this Appendix, we derive the expressions © and © 
for the Laplace transform M(z) of the coherence (0 and 
we show that M(z) is analytic in the upper and lower 
half planes. We first rewrite the expression © as 



M{t) = Tr (e~ iCt p B 



(Bl) 



where £ is a superoperator in the bath Liouville space 
defined by LA = (H B + V)A - A(H B - V) where A 
is any operator in the bath Hilbert space. The Laplace 
transform (7J can then be written as M(z) = Tr[G(z)ps] 
where G(z) — (z — C)^ 1 is the resolvent of the operator 
C. 

We now introduce the superoperators V and Q defined 
by Q = \—V and V A = Tt(A)pb where A is any operator 
in the bath Hilbert space. Since Tr(ps) =1,V and Q are 
projection operators. Using V, Q and (z — C)G(z) = 1, 
we obtain the following coupled equations 



V{z- C)VVG{z)p B 
Q(z-C)VVG(z)p B 



V(z-jC)QQG(z) Pb =Pb 
Q(z - C)Q QG(z)pb =0(B2) 



for the operators VG(z)pb and QG(z)ps- Solving the 
latter for QG(z)pb in terms of VG{z)pb and then sub- 
stituting in the former yields 

V[z-C-CQ(z-QCQ)- 1 QC]VG(z)p B ^PB ■ 

(B3) 

Finally tracing both sides gives [z — "E(z)]M (z) = 1 where 
S(z) is given by 10 • 

We now discuss the analyticity properties of the 
Laplace transform M(z). Consider an eigenoperator A of 
C with the eigenvalue A : CA = \A. The scalar product 
Tt(A^CA) = ATr(AU) can also be written 

Tt(A^CA) = Tt[A^(Hb + V)A-(H b -V)A^aY 

= Tr (A^ CA)* = \*Tr(A Jf A) (B4) 

where the first equality is obtained using the Hermiticity 
of Hb and V and the second one using the cyclic property 
of the trace. Consequently, the eigenvalues of C are real 
and thus the resolvent G{z) and the Laplace transform 
M(z) are analytic in the upper and lower half planes. 

APPENDIX C: SECOND-ORDER SELF-ENERGY 

Here, we show that the second-order self-energy £2 
can be rewritten in terms of the time-dependent correla- 
tion of the interaction operator V as given by 1)10(1. Let 
I a) denote the eigenstates of the Hamiltonian Hb and 
\a,P) = |a)(/3| denote the eigenstates of the correspond- 
ing Liouville operator Lb- Any superoperator T in the 
Liouville space can be expanded in this eigenbasis as 



(CI) 



a, 0,^,5 
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where the scalar product in the Liouville space is defined 
by {A\B) = Tr(A'B). The following decompositions are 
useful for our purpose : 

V = p B £(A 0\ = Y,(a\p B \a)\a, a)(J3, 0\ (C2) 

a, (3 

(z - QCbQ)- 1 = J , p \a,0)(a, (3\ 



a, (3 



Z-E a +Ep' 



£v= [{^\V\l)\a^){ 1 ^\ + { 1 \V\(3)\a,(3){a, 1 \ 



Comparing this expression with the time-dependent cor- 
relation 



(V(t)V) = Tr(p B e ltH Ve- ttH V) 

= J2( a \PB\^)\(a\V\f3)\ 2 e^ E ^ (C4) 

a,0 



where V = 1 — Q and E a is the eigenergy corresponding 
to the eigenstate \a). Using these results we find 



S a (z)= 2^(a|p B |a)|(a|X/|/3)| 2 

a, (3 

1 



(C3) 



z-E, 



Z-Ep+Ea 



we infer that the Laplace transform of 4Re(V(t)V) is 
S 2 (z) as given by flDjl. 
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